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ABSTRACT 



The finned tube contactor of a liquid desiccant (triethylene 
glycol) dehumidif ier system is modeled with a fifth order differential 
equation to simulate the simultaneous heat and mass transfer associated 
with conditioning the air. Potential energy saving aspects of the sys- 
tem are discussed. Computer simulation results presented on psychro- 
metric charts show that significant performance is possible with contac- 
tors of one half the recommended size. The cycle may be regenerated at 
120 to 140°F depending on the system design requirements. 

Figure 4 may be used to estimate system performance. 
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Abbreviations 
LDD 
TEG 

English Letter Symbols 



NOMENCLATURE AND UNITS 

Liquid Desiccant Dehumidification 
Triethylene Glycol 

Cooling area per unit volume for heat transfer, 
ft 2 /ft 3 

Interfacial area per volume for heat transfer. 



Interfacial area per volume for mass transfer. 



c,C 

C,, 



Heat capacity at constant pressure, BTU/lb 
Humid heat, BTU (of water vapor-gas mixture)/ 

(lb dry gas) (°F) defined by Eq. (A. 22) 

Humid heat, BTU (of TEG-water solution) / (lb TEG) 
(®F) defined by Eq. (A. 33) 

Root diameter of tube, ft 

Diffusion coefficient for water vapor in air. 



Mass Transfer coefficient, lb water/(hr) (ft ) 
defined by Eq. (C.l) 

2 

Acceleration due to gravity, ft/hr 
Superficial gas (dry air plus water vapor) mass 



velocity, lb/ (hr) (ft ) 
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G 



max 



G 

s 



h 



G 



h i 

\ 

H 

k 

1 

L 



Gas mass velocity at minimum cross section, 
lb/ (hr) (ft 2 ) 

Superficial solvent (dry air) mass velocity, 
lb/ (hr) (ft 2 ) 

Contact coefficient of heat transfer between gas 
and liquid, BTU/(hr) (ft 2 ) (°F) defined by Eq. (C.2) 
Contact coefficient of heat transfer between gas 
and liquid corrected for simultaneous mass trans- 
fer BTU/ (hr) (ft 2 ) (°F) defined by Eq. (A. 15) and 
(A. 16) 

Contact coefficient of heat transfer between 
liquid and cooling fin, BTU/ (hr) (ft 2 ) (°F) defined 
by Eq. (C.13) 

Conduction coefficient of heat transfer in liquid, 
BTU/ (hr) (ft 2 ) (°F) defined by Eq. (C.ll) 

Enthalpy, BTU/ lb 

Heat of dilution (mixing), BTU/lb 
Thermal conductivity, BTU/ (hr) (ft) (°F) 

Effective fin length, in. 

Superficial liquid (TEG plus water) mass velocity, 
lb/ (hr) (ft 2 ) 

Superficial solvent (TEG) mass velocity, 
lb/ (hr) (ft 2 ) 
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N A 

p * 



Q c 

% 



q 



s 



R 



s 



S 



T 

T 

T 



sh 



sv 

U 



V 



f 



2 

Mass transfer flux of water vapor, lb/ (hr) (ft ) 
defined by Eq. (A. 14) 

Partial pressure in equilibrium with a liquid 
solution, mm Hg 
Total pressure, mm Hg 

Fin pitch, fins/in; Chapter II: pressure 

Heat removed by cooling coil, BTU/hr defined by 
Eq. (A. 28) 

Exothermic heat of dilution, BTU/hr defined by 
Eq. (A. 30) 

Heat transfer flux for sensible heat, BTU/ (hr) 

(ft 2 ) 

Chapter II: universal gas constant; 

Appendix B: bond graph resistance 

Spacing between fins, in. 

2 2 

Interfacial surface, ft /ft cross section 
Thickness of aluminum cooling fin, in; tempera- 
ture, °F 

Chapter II, temperature, °F 
Horizontal tube spacing, in. 

Vertical tube spacing, in. 

Over-all coefficient of heat transfer between the 

2 

liquid and cooling water, BTU/(hr)(ft ) (°F) 

Water velocity through tube, ft/sec 
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x Mass fraction concentration in the bulk liquid, 

lb water/ (lb water + lb TEG) 

X Mass ratio concentration in the bulk liquid, 

lb water/lb TEG 

y Mass fraction concentration in the bulk gas, 

lb water vapor/ (lb water vapor + lb dry air) 

y* Mole fraction concentration in the gas in equilib- 

rium with the bulk liquid, moles water vapor/ 
(moles water vapor + moles dry air) 

y Mole fraction concentration in the bulk gas, 

m 

moles water vapor/ (moles water vapor + moles dry 
air ) 

Y Mass ratio concentration in the bulk gas, lb 

water vapor/ (lb dry air) 

Z Height of finned tube contactor (parallel to 

flow direction), ft. 



Greek Letter Symbols 

T Liquid loading per unit width of fin, lb/(hr)(ft) 

6 TEG solution thickness, ft. 

H Fin efficiency, dimensionless 

X Latent heat of vaporization, BTU/lb 

y Viscosity at bulk temperature, lb/(hr)(ft) 

3 

p Density, lb/ft 

a Ratio of free flow area to face area, dimension- 



less 
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V 


3 

Chapter II: specific volume, ft /lb 


Dimensionless Groups 

Pr 

Sc 


Prandtl number (pc /k) 

P 

Schmidt number (p/p D^) 


Subscripts 

A 


Substance A, water vapor 


AL 


Substance A, liquid water 


B 


Substance B, dry air 


C 


Substance C, TEG 


D 


Substance D, water plus TEG solution; dilution 


f 


Cooling water 


g,G 


Gas 


i 


Interface 


1 ,L 


Liquid 


o 


At the reference temperature 


Bond Graph Symbols 
C 


Capacitance 


C 

m 


Mass capacitance 


°t 


Thermal capacitance 


"a 


Mass transfer flux 


Q 


Heat flux 


R 


Resistance 


R d 


Diffusion resistance 


R 1 


Interfacial resistance to mass transfer 



12 




t 

T 

0 

1 



Reaction resistance 

Heat of dilution flow source 

Latent heat flow source 

Temperature 

Transformer 

Parallel junction, flows sum to zero 
Series junction, efforts sum to zero 
Power bond 

Signal bond 
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I. INTRODUCTION 



Definitions 

Liquid desiccant dehumid i f ica t ion is the reduction of the water 
vapor content of a given volume of air by a liquid sorbent which has the 
property of extracting and holding water vapor brought into contact with 
it. The terms sorbent and desiccant refer to materials with a large ca- 
pacity for moisture compared to their volume and weight. The amount of 
water held by a desiccant will increase or decrease depending on whether 
the vapor pressure of water above the desiccant is less than or greater 
than the partial pressure of water vapor in the air in contact with it. 

Applications 

Dehumidification, whether by liquid desiccants or by other means, 
has important commercial applications. These include preserving ships 
and other surplus equipment which would otherwise deteriorate, and pre- 
venting condensation damage to cargo and internal structures of ships' 
holds which transport dry cargo. 

Liquid desiccant dehumidification (LDD) was investigated in the 
1950 's and found to be operational [7], [22], Currently LDD equipment 
designed for low pressure steam regeneration is commercially available 
and is used in the drying of natural gas prior to pipeline transporta- 
tion. Both hospitals and the food processing industry use LDD not only 
for temperature and humidity control but also for the bactericidal prop- 
erties of the family of glycol desiccants. Accepted high equipment cost 
prevents wider application of LDD for the general purpose of space de- 
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humidification or space cooling. These costs are on the order of $10 
per square foot of contactor in typical LDD systems which remove about 
500 pounds of water per hour with dehumidification and regeneration con- 
tactor areas of 5000 sq. ft. each [27]. 

Potential Savings 

If the potential to regenerate an LDD cycle at low temperatures is 
exploited and waste heat or solar energy is utilized, then potential 
energy savings can be applied against acquisition cost. Such cost may 
be minimized by rational rather than empirical equipment selection. 

A review of current space dehumidification and cooling practice 
will indicate the potential energy savings of LDD. 

In conventional air conditioning systems the air is cooled signifi- 
cantly below the desired temperature to effect moisture removal. (Refer 
to Figure 1). 



CHILLED WATER ENERGY 



i L 



WARM. MOIST AIR 


COOLING 


VERY COOL, DRY AIR 


REHEAT 


COOL DRY AIR 




COIL 




COIL 





T 



Figure 1. Conventional Air Conditioner with Reheat 
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Additional energy is used as this very cool, dry air is then reheated 
to the desired space temperature. The popularity of this system stems 
from the previously low cost of energy and the simple control system re- 
quired. One dependent variable, outlet humidity ratio, is uniquely con- 
trolled by the chilled water temperature through the condensation of ex- 
cess moisture above the dew point temperature. The other dependent vari- 
able, desired space temperature, is now uniquely controlled at a constant 
humidity by the quantity of reheat energy added. 

The control system complexity increases when an attempt is made to 
save reheat energy by mixing the very cool dry air with warm moist air 
bypassed around the cooling coil. (Refer to Figure 2.) 



CHILLED WATER 




Figure 2. Conventional Air Conditioner with Bypass Mixing 

Both of the two independent variables, chilled water temperature and bypass 
air flow rate, simultaneously modify both of the dependent variables, 
outlet air temperature and humidity ratio. 

The LDD system can save the energy of excess chilling for moisture 
removal through its ability to dehumidify at ambient temperatures. Re- 
moval of the latent heat of condensation is the primary function of the 
cooling water. (Refer to Figure 3.) 



16 
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WARM MOIST AIR AIR-LIQUID COOL DRY AIR 



CONTACTOR 




DILUTE 

LIQUID DESICCANT 



Figure 3. Liquid Desiccant Dehumidification System 



The dehumidif ication contactor of the LDD system, a three fluid heat ex- 
changer in which mass transfer between two of the fluids occurs, can 
supply cool, dry air without further heating or cooling. Thus, the re- 
heat energy can also be saved. 

In the LDD contactor the two dependent variables, air temperature 
and humidity ratio, are controlled by the independent variables, liquid 
desiccant concentration and cooling water temperature. A tradeoff for 
the potential energy savings of excess chilling and unnecessary reheat 
is a control system of sophistication similar to that of the bypass sys- 
tem shown in Figure 2. 

The LDD cycle compares favorably to an alternative technique for 
space cooling such as the lithium bromide absorption cycle. While the 
absorption cycle is closed and the desiccant cycle is open, both systems 
have coefficients of performance on the order of 0.6 [11]. The absorp- 
tion cycle requires regeneration temperatures around 200°F, while the 
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LDD cycle can be regenerated at temperatures around 140 °F. Of these 
two systems the absorption system has been the subject of considerable 
analysis and simulations [9]; however, there is no well-established tech- 
nique for space cooling or dehumidification with low temperature regener- 
ation. 

Purpose 

In view of the above, a model for analysis of the liquid desiccant 
dehumidifier contactor of Figure 3 is presented as the basis for a ra- 
tional design. The primary function of the model is to assure that the 
contactors are not over-specified so that acquisition costs may be kept 
as low as possible. The contactors (dehumidification and regeneration) 
represent a significant portion of the total cost of an LDD system. 

Empirical equipment selection of the contactors is standard prac- 
tice and can result in overspecification to ensure proper operation. 
Manufacturers' application recommendations specify equipment which oper- 
ates with the exit conditions of air and liquid temperatures, humidity 
ratio, and desiccant concentration at equilibrium. 

By contrast, rational design implies more than black box modeling 
of the contractor. The computer model of the contactor which was devel- 
oped is based on a complete set of energy balances, material balances, 
rate equations, and equilibrium relationships. Through computer simu- 
lation of the heat and mass transfer in the contactor, the rate at which 
the temperatures, humidity ratio, and desiccant concentration approach 
equilibrium as a function of flow distance can be studied. With a speci- 
fic design the contactor flow length and associated cost can be compared 



18 



to the incremental heat and mass transfer of additional contactor flow 
length increments. 
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IX . PSYCHROMETRY 



The Psy chrome trie Chart 

Here psychrometry is concerned with the determination of the prop- 
erties of the mixture of air and water vapor as the mixture flows through 
the contactor. These properties are conveniently diagrammed in an abbre- 
viated psychrometric chart. Figure 4, which plots dry bulb temperature, 
humidity ratio, wet bulb temperature, saturation temperature, and rela- 
tive humidity for the air -water vapor mixture. 

Dry bulb temperature represents the air and liquid contact tempera- 
ture and is the abscissa. 

Humidity ratio (H.R.) is the pounds of water carried by one pound 
of dry air. Alternatively, grains per pound of dry air (grains) is used 
where 7000 grains equal one pound of water. The ordinate shows both 
humidity ratio and grains. 

Relative humidity is the ratio of partial pressure of the water 
vapor to the vapor pressure of water at the air temperature and is indi- 
cated by the dotted lines. Relative humidity is expressed in percent. 

Wet bulb temperature lines run upward and to the left to the satu- 
ration curve and very nearly represent the path followed by adiabatic 
(spray) cooling. 

Dew point temperature is the temperature to which the mixture must 
be cooled to become saturated and is found by following a horizontal line 
to the saturation temperature curve. 

Saturation temperature represents the maximum moisture the air can 
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Figure 4. Triethylene Glycol Equilibrium Humectant Values on Psychrometric Chart 
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hold before fogging ensues. 

Enthalpy and specific volumes are plotted on more detailed charts 
[26], Any two of these parameters specify the air-water vapor mixture 
and allow all the others to be determined. The remaining parameter, 
humectant value, is discussed later under its own section. 

Dry bulb temperature and humidity ratio were chosen as the two de- 
pendent variables of interest because: (1) dry bulb temperature is 

easily measured and has significance to people; and (2) the humidity 
ratio directly represents the water content of the air and is based on 
one pound of dry air which is a constant in any psychrometric process. 

Psychrometric Chart Processes 

The process of heating or cooling moist air above the dew point is 
represented by a horizontal line since the humidity ratio is unchanged. 
This process is known as sensible heating or cooling. 

The process of cooling below the dew point entails condensation of 
vapor and is represented by a horizontal line extending left to the satu- 
ration curve where the process follows the curve to a lower humidity 
ratio. This is roughly the dehumidif ication process of conventional air 
conditioning. 

A downward vertical line on the chart represents an isothermal dehu- 
midification process or latent cooling. This may be accomplished directly 
by the LDD system with cooling to remove only the latent heat of conden- 
sation. 



Liquid Desiccants 

A liquid desiccant removes water vapor from the air at a rate direc- 
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tly proportional to the difference between the partial pressure of the 
water vapor immediately above the liquid desiccant. The partial pres- 
sure of water vapor in the air is reduced and the liquid desiccant is 
diluted in the process until equilibrium occurs. The process may be re- 
versed by heating the desiccant. That is, the liquid desiccant is con- 
centrated as the partial pressure of water vapor in the air increases. 

The vapor pressure of the liquid desiccant is a function of temperature 
and concentration. 

Humectant Values 

Equilibrium humectant values are the percent by weight of a liquid 
desiccant in water solution required to maintain vapor pressure equili- 
brium in contact with air of various temperatures and humidity ratios. 

The humectant values for triethylene glycol (TEG) were added to the 
standard psychrometric chart to produce a convenient diagram revealing 
the relationship between TEG concentration and the air-water vapor mix- 
ture at equilibrium. The humectant values on Figure 4 were developed 
from graphs of the vapor pressure of water over triethylene glycol of 
various weight concentrations and temperatures [28] [37]. These equilib- 
rium humectant values give an indication of the TEG concentration necess- 
ary to dehumidify to a specified outlet condition with a constant TEG 
concentration. 

The solid lines which are the TEG equilibrium humectant values rep- 
resent slightly different TEG weight concentrations in the 90% to 99% 
range depending on the data base: Union Carbide [37] or Dow Chemical [281 

A third source [12], presented data for only three concentrations, 40%, 
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60%, and 80% and these were consistent with the first two sources within 
5%. The vapor pressures presented by Union Carbide and Dow Chemical 
differ by approximately 20% at 90% concentration and 100% at 98% concen- 
tration. Union Carbide data presented the higher vapor pressures. The 
more conservative humectant values during dehumidification are derived 
from the Union Carbide data while the more conservative humectant values 
during regeneration are derived from the Dow Chemical data. 

In the area of Figure 4 in which dehumidification takes place, 70 
to 80 °F and 0.010 H.R., the change in TEG humectant value with respect 
to temperature is about 1%/°F. At a 0.025 H.R. and 120 to 130 °F where 
regeneration can take place, the ratio is about 0.2% /°F which indicates 
that TEG regeneration should not be particularly sensitive to regener- 
ation temperatures. 

Vapor Pressure Prediction 

An equation for the equilibrium partial pressure above the liquid 
desiccant is necessary for the computer simulation. In addition to 
accuracy, a suitable equation will have continuous derivatives over its 
range and will be relatively efficient with respect to calculation time. 
Two relationships with these characteristics which predict vapor pres- 
sure as a function of TEG contact temperature and TEG weight concen- 
tration were developed with the aid of a curve fitting program and a 
programmable calculator: 

From Dow Chemical data. 



p* [mm Hg] = exp 



-9346 



+ 20.5749 - 




21.36(98-P X ) 98-P t 



T + 460 




( 2 . 1 ) 



0.28294 
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From Union Carbide data. 



p* [mm Hg] = exp 



-9259 



+ 20.407 - 1.4333(: 



—1 ) 1‘5 + 

100.5-P t ; 



T + 460 




( 2 . 2 ) 



0.26921 



where T = contact temperature, °F 



P, 



T = TEG weight fraction concentration, %. 



Both relationships are accurate within about 3% within the ranges 
W < T < 300, 75 £ C £ 99 and are slightly less accurate for 0 < C <; 75, 
Eq. (2.2) was used exclusively in the simulation. 

Vapor Pressure Correlation 

The method used in synthesizing equations (2.1) and (2.2) is based 
on the Claus ius-Clapeyron equation which relates the slope of the vapor- 
pressure curve to the latent heat of vaporization 



dP _ A 



(2.3) 



dT T(V g - V L ) 



where = specific volume of the saturated vapor 
= specific volume of the saturated liquid 
A = latent heat of vaporization 



Using the ideal-gas law and neglecting in comparison with 




(2.4) 



If A is fairly constant over the temperature range 



In P =- — + constant 



(2.5) 
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Equation (2.5) suggests that a plot of log P against ^ will be linear 
over small temperature ranges. Antoine's equation published in 1888 
allows for slight variation of X with temperature [26] 

lnP = a ' t~+T (2 - 6 > 

and is probably the best vapor-pressure correlation equation for a simple, 
three-constant (a, b, and c) equation. Equation (2.6) is the basis for 
equations (2.1) and (2.2). 
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III. LIQUID DESICCANT DEHUMIDIFICATION SYSTEM 
Operating Cycle 

One system for liquid desiccant dehumidification (LDD) is shown 
in Figure 5. The dehumidifier section on the left accepts warm moist 
air from the atmosphere and contacts the air with the liquid desiccant 
in the dehuraidifier contactor. Intimate contact between the air and 
desiccant over a large area occurs through the formation of a thin film 
of desiccant as it is sprayed over the extended surface contactor. The 
heat of condensation and the heat of dilution flow through the liquid to 
the cooling water which can be sea water, well water, or water from a 
cooling tower. The quantity of cooling water required is strongly de- 
pendent on the sensible and latent heat removed from the air by the 
liquid desiccant and weakly dependent on the desiccant circulation rate 
between the dehumidifier and regenerator. 

The air leaves the contactor flowing downward with entrained desic- 
cant and loses most of the liquid to the sump while crossing it. The re- 
mainder of the entrained liquid is eliminated as the cool, dry air flows 
upward to the controlled space. 

The dry bulb temperature of this cool, dry air is weakly dependent 
on the water removed from the air but is strongly dependent on the temp- 
erature of the liquid desiccant. Accordingly, the air flowing through 
the dehumidifier can be warmed or cooled as it is dried depending on the 
temperature of the liquid desiccant and the temperature of the cooling 
water. 



Concentration of the liquid desiccant is more easily achieved with 



COOL, DRY HOT. WET 

AIR TO SPACE AIR OUT 
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Figure 5. Four Pump Liquid Desiccant Dehumidification System 
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dry air which would normally be exhausted from the dehumidified space, 
but can also be achieved with outside air. The heat of vaporization 
and dilution is supplied through the heating coil in the regenerator 
contactor to the liquid desiccant as it contacts the relatively dry air. 
Heating evaporates the liquid desiccant though at a much reduced rate 
than that of the water. A reflux coil maintained at about 90 °F is used 
as a condenser for the desiccant to prevent excessive loss in the hot 
exhaust. 

Figure 5 also shows the four pumps necessary for the liquid desic- 
cant movement — two spray pumps and two feed pumps. The liquid-to-liquid 
heat exchanger between the conditioner sump which is slightly warmer 
than the cooling water, and the concentrator sump which is slightly 
cooler than the heating water, improves the system thermal efficiency. 

An air-to-air heat exchanger between the hot wet exhaust and the cooler 
inlet air of the regenerator will also increase the system thermal effi- 
ciency. 

Figure 6 shows a similar system with a combined sump which requires 
only one pump. While this system is more compact and uses fewer pumps, 
thermal efficiency can be degraded since the sump operates at a mean 
temperature. With equa.1 dehumidifier and regenerator desiccant recircu- 
lation rates the sump operates as a counter flow heat exchanger with an 
effectiveness of 0.5, or less than is normally expected from the counter 
flow heat exchanger of Figure 5. 

In either system the solution in the sump will require treatment as 



it will concentrate solids absorbed from the air. 



COOL. DRY HOT. WET 
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Figure 6. One Pump Liquid Desiccant Dehumidification System 
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Air-Desiccant Contactor 

The contactor is the part of the system in which the air is brought 
in contact with the liquid desiccant for dehumidification or regener- 
ation. The contactor should have high heat and mass transfer coeffici- 
ents but yet have relatively low air flow resistance. The contactor 
should also be light, compact, inexpensive and retain its heat and mass 
transfer characteristics after extended use. 

Although the contactor may consist of a simple spray with the re- 
sulting air-desiccant interface being the surface of droplets, entrain- 
ment would be excessive since the finely divided liquid is difficult to 
separate from the air. In addition, it would be difficult to remove the 
heat of condensation and dilution. A drop spray onto extended surface 
finned tubes. Figures 7 and 8, prevents atomization of the liquid desic- 
cant and has a better mass transfer capacity per volume than a simple 
spray [33]. 




Figure 7. Four Row Finned Tube Contactor Dimensions 
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Dimensions for Figure 7 

T = 2.00 in T = 1,73 in 

sh sv 

Width = 48 in Length = 127 in 

Height = 6.69 in/ (four rows) 




P = 5 fins per inch = 0.75 in 

t = 0.010 in T , = 2.00 in 

sh 



Figure 8. Dimensions for Determining Free Flow Area of Contactor 
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The cooling water flows through the tubes in a parallel-counter- 
flow configuration with an even number of tube passes, while the air and 
liquid desiccant flow over the extended surface. Thus, the dehumidifi- 
cation contact surface and cooling coil are conveniently packaged in one 
unit. The cooling water serves to: (1) keep the desiccant cool and thus 

lower the vapor pressure; and (2) cool the air. 

The regenerator contactor has the same physical characteristics 
and performance criteria as the dehumidifier contactor, although the size 
may differ because of different air flow and liquid desiccant recircu- 
lation rates. Therefore, the computer model presented later simulates 
the heat and mass transfer of both the dehumidif ier and regenerator con- 
tactors. 

The air flow can be either parallel with or counter to the liquid 
desiccant flow direction, depending on the particular application. Due 
to the dilution of the liquid desiccant as it flows through the contac- 
tor, the exit humidity ratio is not as low as can be achieved with the 
counter flow arrangement. However, the downward parallel flow configu- 
ration does not inhibit spraying the liquid on the contactor and allows 
some separation of the liquid and air by a reversal of flow direction. 

The parallel flow configuration is best suited for the processing of 
large volumes of air, while the counter flow configuration can deliver 
the very dry air needed in manufacturing operations such as the handling 
of hygroscopic materials. The parallel flow configuration represents 
our physical system. 
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Other System Equipment 

Other equipment items, such as the main circulating blowers or 
liquid sorbent pumps, should be selected on the basis of good engineer- 
ing design with cognizance of the difficulty of containing TEG because 
of its propensity to leak past defective fittings or packings. 

I.DD System Advantages 

1. Low temperature regeneration as opposed to other thermal cycle air 
conditioners is possible and allows more efficient operation of solar 
collectors. 

2. Regeneration is not very sensitive to heating water temperature, un- 
like other thermal cycle air conditioners. 

3. Solar or recovered waste power can be used to regenerate the liquid 
desiccant . 

4. The system is simple in construction and operation since the cycle 
is without high pressure liquids or gases. 

5. Precise control of air temperature and humidity results from a speci- 
fic TEG concentration and cooling water temperature. 

6. Co-sorption benefits include the removal of bacteria, odors and 
various airborne irritants [1], 

7. The expense of a separate cooling coil can be saved or a smaller 
after cooler can be specified if chilled water is used in the dehumidi- 
fication contactor. 

8. Excess dehumidification capacity when obtained from low cost waste 
heat or solar energy may be advantageously converted to sensible cooling 
by spray cooling. The spray process cools and humidifies essentially 
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along the lines of constant wet bulb temperature in Figure 4. 

9. In ocean engineering applications dew point temperatures on the 
order of 10 °F below the structure temperature are required to prevent 
condensation [10]. The combination of LDD and cold sea water can produce 
sub-freezing dew point temperatures as the antifreeze qualities of TEG 
prevent contactor frost accumulation. 
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IV. METHOD OF ANALYSIS 

Mathematical Model 

The analysis of the contactor is based on the development of a 
mathematical model suitable for a computer simulation. The substance 
of the model developed is presented here; detailed information may be 
found in the appendices. 

The model of the contactor is a fifth order differential equation 
composed of three first order differential equations for heat transfer 
and two first order differential equations for mass transfer. The 
three heat transfer state variables are the gas, liquid, and cooling 
water temperatures as shown in Figure 9. 
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Figure 9. Locations of Representative Temperatures and Sources 



36 



TEG Solution thickness 0.007 to 0.010 in 



Fin thickness 



0.010 in 



Tube wall thickness 



0.035 in 



Tube diameter 



0.750 in O.D, 



Dimensions for Figure 9 



The two mass transfer state variables are the masses of the gas, (water 
vapor plus dry air), and liquid, (water plus TEG). 

All differentials are taken with respect to flow distance in a 
steady state model. This is a reasonable assumption, as start up trans- 
ients will be negligible within a minute or two of operation from a de- 
termination of the system RC time constants. 

Actual variables used in the program for the mass state variables 
are the superficial mass velocities of the gas and liquid. Thus, the 

five independent variables in this model are t , t, , t,., G, and L. All 

8 1 f 

the other variables are dependent variables and can be calculated from 
known relationships with knowledge of the five independent variables. 

Two distinct methods were used to develop the mathematical model. 
Appendix A presents a derivation using the traditional method for analyz- 
ing mass transfer operations based on material and energy balances in 
control volumes of a differential section of the contactor. Appendix B 
uses a general method. It develops a bond graph model of the contactor 
from which the governing differential equations rigorously follow. The 
equations developed by the two methods are equivalent. 
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The five diffenential equations from Appendix A and B are: 



dt 

_ a = 

dZ 



R C 
8 g 



(4.1) 



dZ 
dt _ 



f - Vm 



dL , 
dZ 



Vm 



N (1094 + 0.445t ) + + — 

A S K g K 1 



1 8 



(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 



Simulation 

These five differential equations and the other supporting equa- 
tions listed in Appendix B which calculate variables such as the mass 
transfer flux, humidity ratio, and TEG weight fraction concentration 
form a subroutine. A master program uses the forth order Runge-Kutta 
method and this subroutine to produce an incremental numerical solution. 

Inputs number 19. The geometry of the contactor is specified with 
eleven inputs. Initial temperatures take three inputs. Flow rates and 
concentrations each take two inputs and a final input controls the in- 
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crement size and the stability of the solution. 

A substantial portion of the subroutine calculates initial values 
of approximately 100 variables including the fluid properties and trans- 
fer parameters discussed in Appendix C. 

A minor portion of the subrountine is repeatedly called during the 
actual solution. Even so, approximately three minutes are required to 
calculate the 2500 increments of a 16 row contactor simulation. Each of 
the five state variables is calculated four times per increment while 
the auxiliary variables are calculated only once per increment. 

The primary outputs of interest are the state variable t^ and the 
auxiliary variables: Y, humidity ratio; P^, percent TEG weight fraction 

concentration; and Z, the distance in air flow direction. Cross plots 
of t^ and Y produce operating lines on psychrometric charts while cross 
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V. RESULTS OF SIMULATION 



Dehumidifier Contactor 

In assessing the performance of the dehumidifier contactor, the 
operating line on a psychrometric chart furnishes the pertinent infor- 
mation: the inlet and exit conditions and the path between them. The 

psychrometric chart. Figure 10a, displays a parametric study of dehumidi- 
fier operating lines for three TEG inlet concentrations, 92%, 94% and 
96% and for three cooling water temperatures, 60°, 70° and 80°F. All 
simulations begin at 92 °F dry bulb and 76 °F wet bulb (approximately 
50% relative humidity) which is a typical outdoor summer temperature 
specification for performance rating of air conditioning units. 

Marks placed on the operating lines of Figure 10a represent equal 
flow distances of four rows (6.69 inches) except for the special mark 
which represents conditions two rows (3.34 inches) from the inlet (top 
of figure). The parameters specified in the figure are the manufacturer’s 
recommendations for a 16 row contactor. For all nine simulations the 
mark indicating 8 rows is very close to the 16 row mark (within 2 to 4°F 
and 3 to 4 grains) , and represents a 50% savings in contactor acquisition 
cost. Also, a 4 row contactor using 60°F cooling water is nearly equiva- 
lent in performance to a 12 row contactor using 70°F cooling water for 
the same TEG inlet concentration. Such equivalent performance points 
can be the basis for optimizing cooling water operating costs and equip- 
ment acquisition costs for a specific application. 



The dilution vs. flow distance of the TEG solution for the Simula- 
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Figure 10a. Dehumidif ier Operating Lines for Three TEG Concentrations 

and Three Cooling Water Temperatures on Psychrometric Chart 
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tions of Figure 10a is presented in Figure 10b. In Figure IQb, 
for the same inlet TEG concentrations, the lowest exit concentrations 
correpsond to the lowest temperature cooling water simulations and the 
drier conditioned air. The data on Figures 10a and 10b are from the 
same simulations. 

A correlation of Figures 10a and 10b with the equilibrium humecta,nt 
values of Figure 4 shows that the exit conditions of the 16 row contac- 
tor, operating as recommended by the manufacturer, are at equilibrium. 

For example, from Figure 10b an exit TEG solution concentration of 90% 
corresponds to an inlet concentration of 94% with 60°F cooling water. 

On Figure 10a the 94% simulation with 60°F cooling water ha,s exit con- 
ditions of 60.5°F and 0.005 H.R. or 35 grains. On Figure 4 a dry bulb 
temperature of 60.5°F and a 0.005 H.R. corresponds to a TEG equilibrium 
humectant value of 90% (the Union Carbide data base is used in the compu- 
ter simulation), which is the exit TEG concentration of Figure 10b. 

When psychrometric processes are graphically displayed in the liter- 
ature the path lines are often straight. The path lines for the liquid 
desiccant dehumidification simulation are not straight and can be divided 
into two regions (Refer to Figure 10a.) The first region is nearly all 
latent cooling (isothermal dehumidification) in which the ability of the 
cooling water to perform sensible cooling (dry bulb temperature change) 
is inundated by the exothermic condensation process. In the second reg- 
ion sensible cooling is more pronounced. 

The first region physically represents the flow through approxi- 
mately the first two rows of a 16 row finned tube contactor or approxi- 



DEHUMIDIFIER 48" X 127" 
G WATER TEMPERATURES: 
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Figure 10b. Dehumidifier TEG Concentration vs. Flow Distance 
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mately 3.34 inches through a contactor 26.75 inches high. The combin- 
ation of a high rate of mass transfer at the inlet (see Figure 11) and 
the associated heat of condensation partially explain the nearly iso- 
thermal region. The exothermic heating due to mass transfer is obvious 
in Figure 12 which shows the variation of the liquid, interface, and air 
temperatures against distance. 

Chemical engineers frequently use a graph of the operating line and 
equilibrium curve in analyzing the absorption of a substance in gas- 
liquid operations. Figure 13 shows these lines for the present system. 
Some simplifications of the calculation procedure can be made if the 
two lines are straight or are straight and parallel [23], However, for 
the LDD system these lines are not straight because of the temperature 
influence on the equilibrium vapor pressure and hence a computer simu- 
lation is required. 

The operating line of Figure 13 represents the local mole fraction 
of water vapor in the air at the local mole fraction of water in the TEG 
solution. The equilibrium curve represents the equilibrium mole fraction 
water vapor in air at the local mole fraction water in the TEG solution. 

The simulations for Figure 10a were iterated to achieve identical 
inlet and exit TEG solution temperatures. This corresponds to recircu- 
lating the solution through the contactor. Thus, the sump temperature 
for each simulation is slightly less than the air outlet temperature. 

The greatest change in liquid temperature is near the inlet (refer to 
Figure 12) so that the liquid temperature at the 8 row point of a 16 row 
contactor differs little from the liquid temperature at the exit. For 
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Figure 13. Operating Line and Equilibrium Curve for Dehumidifier Contactor 
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this reason the 8 and 12 row points may be accepted as exit conditions 
even though they are not. Figure 14 verifies this point as the 8 row 
exit marks obtained by iterating the simulation for equal inlet and exit 
liquid temperatures nearly coincide with the 8 row marks of Figure 10a, 
the 16 row contactor. 

Figure 15 shows the results of an air flow rate which is 75% of the 
air flow rate of the Figure 14 simulations. Slightly drier, 1.5 grains, 
and slightly cooler, 1°F, air is the result of a reduced air flow rate 
on the exit conditions. 

Specifying a higher aspect ratio (defined as a greater face area to 
flow distance ratio for a given volume of contactor) will allow the 
treatment of more air. Figures 10a and 14 support this statement. For 
example, the 16 row contactor with a 48" by 127" face area could be 
arranged as an 8 row contactor with a 96" by 127" face area. This high- 
er aspect ratio contactor will process twice as much air assuming con- 
stant TEG recirculation rates and air flow rates per unit face area to 
exit conditions within 2 to 4 °F and within 2 to 4 grains of the 16 row 
contactor exit conditions. 

An indicative result of an 8 row dehumidifier contactor (from Fig- 
ures 4 and 10a) is that with 70°F cooling water and 92% inlet TEG con- 
centration the inlet air at 92°F and 111 grains exits at 72.5°F and 56 
grains or a dew point temperature of 51°F. 

Regenerator Contactor 

In assessing the performance of the regenerator contactor the exit 
TEG concentration is the significant parameter. Figure 16a shows the 
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Figure 15. Dehumidifier Operating Lines at 75% Air Flow Rate 
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ability of the regenerator contactor to concentrate 90% and 92% TEG so- 
lution at temperatures of 120°F and 140°F. At constant temperature re- 
generation the exit concentration variations are less than the inlet 
concentration variations. For example, at 140°F regeneration, the in- 
let concentration variation for two simulations is 2% while the exit 
variation is 1/2%. Figure 16b shows the regenerator operating lines on 
a psychrometric chart. Both Figures 16a and 16b are for the same simu- 
lations. 

The regenerator exit conditions are nearly at equilibrium as a 
study of Figures 16a, 16b and 4 will verify. For example, the nearly 
95% exit concentration of Figure 10a for 90% TEG inlet concentration and 
120°F heating water corresponds in Figure 10b to 118°F and 0.021 H.R. or 
147 grains. On Figure 4 a dry bulb temperature of 118°F and 0.021 H.R. 
corresponds to an equilibrium humectant value of about 95% (Union Carbide 
data base). The second equilibrium relationship, equation (2.2), for a 
contact temperature of 118*^ and 95% TEG concentration predicts a vapor 
pressure of 24.6 mm Hg which corresponds to 0.0208 H.R. 

Y = 0.622 (—2 ) = 0.622(^14) = °* 208 (6.1) 

/bU - p /O0.4 

The significant result of the 8 row regenerator contactor simulation 
is that TEG concentrations of 97% to 97.5% can be achieved with 140°F 
heating water, and 95% to 96% TEG concentrations can be achieved with 
120°F heating water. 
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Figure 16a. Regenerator TEG Concentration vs. Flow Distance for Two TEG 
Inlet Concentrations and Two Heating Water Temperatures 
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Figure 16b. Regenerator Operating Lines on Psychrometric Chart 
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VI. CONCLUSIONS 

1. Contactor equipment costs may be substantially reduced because half- 
sized contactors can provide 95% of the dehumidification and 90% of the 
regenerated TEG concentration as compared to the full sized contactors 
recommended by the equipment manufacturers. 

2. Most effective dehumidification per unit contactor volume occurs at 
the highest ratio of face area to flow distance. 

3. The sensitivity of the regeneration process to the regeneration 
temperature is minimal. 

4. Performance estimates of dehumidification processes which continue 
to equilibrium can be prepared from Figure 4, Triethylene Glycol Equili- 
brium Humectant Values on a Psychroraetric Chart. First, assume that the 
exit conditions lie on a vertical line at the cooling water temperature. 
Then, since each point on the line represents a humidity ratio and a 
humectant value, and the sum of the air and liquid masses is constant, 
the equilibrium point can be found. 
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APPENDIX A 

CONTROL VOLUME DERIVATION OF RATE RELATIONSHIPS 

The equations derived in this appendix are concerned with the inter- 
phase transfer of mass (water vapor) and energy which occurs during de- 
humidification of the humid air by the concentrated liquid desiccant. 
Knowledge of the equilibrium characteristics of the system is necessary 
for an understanding of the mass transfer operation and, because of 
simultaneous heat transfer, enthalpy characteristics must also be con- 
sidered . 

Operations such as liquid desiccant dehumidification quickly become 
complex owing to the mass transfer in both the gas and liquid phases and 
the large heat effects which accompany condensation and dilution. Fur- 
thermore, the need to remove the heat of condensation results in a non- 
adiabatic operation. 

Figure 17 represents the direct contact of the gas (dry air and 
water vapor) and liquid (TEG and water) flowing side by side in parallel 
flow and separated by the gas-liquid interface. The changes in tempera- 
ture, concentration, etc., are all differential over the section. 

Based on Figure 17, the balance of this appendix develops: 

(1) An overall material balance for control volume III; 

(2) Energy balances for control volumes I, II and III; 

(3) An interface control volume energy balance. 

Next, parameters are separated into catagories and finally, a solution 
sequence is presented . 
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Figure 17. Differential Section of Contactor 
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The interfacial surface area in the section is dS. For the finned 
tube contactor the heat transfer area per unit volume, a^, is known 
from calculations. Because of the low surface tension of the TEG solu- 
tion complete wetting of the heat transfer area may be assumed so that 
the mass transfer area per unit volume, a^, is equal to the heat trans- 
fer area: 



*H = “M (A - 1) 

and, 

dS = a M dZ (A. 2) 

where the volume per unit section is dZ. 

Consider the mass transfer operation of Figure 17 as conducted in 
a steady state cocurrent fashion in which Lg and Gg represent the super- 
ficial mass flow rate of the non-diffusing pure TEG and dry air respec- 
tively. Since the values Lg and Gg are constant through the dehumidifier, 
they are the basis for the material balance. 

At the entrance to the dehumidifier the gas phase contains G total 
lb/(hr)(sq ft) consisting of non-diffusing G g lb/(hr)(sq ft) and the 
diffusing water vapor. The water vapor concentration is given by Y the 
mass ratio, or equivalently y the mass fraction. The relationships 
among these variables are: 

Gg = G(l-y) = y VY~ (A ’ 3) 




(A. 4) 
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y = 



Y 

1 + Y 



(A. 5) 



Similarly, the liquid consists of L total lb /(hr) (sq ft) containing 
x mass fraction of water (or mass ratio X) and Lg lb/(hr)(sq ft) of non- 
diffusing TEG. The relationships among these variables are: 



L s - L(1 - x) ■ T7I 



(A. 6) 



X = 



x 

1 - X 



(A. 7) 



x 



X 

1 + X 



(A. 8) 



Control Volume III Material Balance 

The quantity of water vapor in the gas is Gy lb/(hr)(sq ft). In 
the liquid it is Lx lb/(hr)(sq ft). Hence the overall material balance 
about control volume III of Figure 17 for the water is: 



In = Out 

Lx + Gy = (L+dL)(x+dx) + (G+dG) (y+dy) (A. 9) 

expanding, 

Lx + Gy = Lx + Ldx + dLx + dLdx + Gy + Gdy + dGy + dGdy (A. 10) 

neglecting the second order terms 

d(Lx) = -d(Gy) (A. 11) 

or 



LgdX = -G s dY 



(A. 12) 
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The mass transfer rate in lb vapor/(hr) (sq ft) of dehumidifier cross 
section) is [36] 

Wi ' 'V Y (A - 13) 



where 



lb water vapor absorbed 
(hr) (interfacial surface) 



F G ln( 




-) 



(A. 14) 



Here y* is the equilibrium mole fraction concentration of water vapor 
over the TEG solution at the interface; y is the bulk mole fraction 
concentration of water vapor in the gas and is the mass transfer co- 
efficient for large mass transfer rates and is defined by Eq. (C.l). 
Equation (A. 14) is a special case for the absorption of one component. 

The gas sensible heat transfer rate across the interface in BTU/ (hr) 
(sq ft) is [36] 

^sgV 2 " l -e^-V - /h G T (t G - C i> V 2 (A ’ 15 > 

or 

q sG a H dZ = h G (t G ' W 2 (A ‘ 16) 

where = heat capacity of water vapor at constant pressure, 

BTU/ (lb mass) 

hg = convection heat transfer coefficient, BTU/(hr)(°F) 

(sq ft interface) and is defined by Eq. (C.2). 

Here h^ includes the effect of mass transfer on heat transfer by 
accounting for the sensible heat carried by the diffusing water vapor. 
Radiation has been neglected. 
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The sensible heat transfer through the liquid in BTU/(hr)(sq ft) 
is [36] 

VV Z - h L a H <t i - C L )dZ tA ‘ 17) 

Energy balances around each of the three envelopes of Figure 17 can 
now be developed. 



Envelope I Energy Balance 

In = Out 

G S H G - G S <H G + dH G> + %G% iZ + W Z(C i ( VW (A - 18) 



where the last term is the enthalpy of the transferred vapor relative 

to a reference state at t Q where is the latent heat of vaporization 

at t . Since 
o 



N A a M dZ = -d (Gy) 




(A. 19) 



substitution yields 

0 ’ G S dH G + h G a H Ct G- t i )dZ - G S (C A <t G- t o )+ V dY < A ' 20 > 



The enthalpy of the air-water vapor mixture per pound of dry air is 
H G * C B <t G' t o > + Y<C A < V t c )+X o ) <A - 21) 



where C R = heat capacity of dry air at constant pressure, 

BTU/(lb)(°F) 

Define 

C S - C B + YC A <A - 22) 

= humid heat of vapor-air mixture, BTU/(lb dry air) (°F) 
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so that 



Differentiating equation (A. 21) 



dH G - C B dt G + YC A dt G + ^ C A^ C G _t o^ + V dY 

or, 

dH G - C S dt G + (C A ( V t o ) + V dY 

Combining with (A. 20) gives the envelope I energy balance 

- G S C S dt G ' VH W “ Z 
Envelope II Energy Balance 

In + heat generation = Out 

l sVV\V ZC AL (t r t o ) + W Z ‘ VW + Q C 

Here 

Q C * U W t f )dZ 

= heat removed by cooling coil, BTU/hr 
where U = over-all heat transfer coefficient between the 

liquid and cooling water, BTU/(hr)(sq ft)(°F) 
a^ = cooling heat transfer area per unit volume 
sq ft/cu ft 

t f = cooling water temperature, °F 



(A. 23) 

(A. 24) 
(A. 25) 

(A. 26) 

(A. 27) 
(A. 28) 



Substituting equations (A. 17) and (A. 19) into (A. 27) yields 
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% - G S C AL (t i- t o )dY + WW" = L S dH L + Q C (A ' 29) 

Here 

q d = W zh d = - G sV Y (A - 3 °) 

= the exothermic heat of dilution (mixing) , BTU/hr 
where = heat of dilution, BTU/(lb water transferred) 

Substituting equations (A. 28) and (A. 30) into (A. 29) yields 

- G S (H D +C AL (t i- t o» dY + h L a H (t r t L )dZ * W Ua C « V t f )dZ < A ' 31 > 

The enthalpy of the TEG - water solution per pound of TEG is 

«L ' WV + XC AL (t L- t o> (A ' 32) 

where C c = heat capacity of TEG at constant pressure 
BTU/(lb)(°F) 

= heat capacity of liquid water at constant pressure 
BTU/(lb)(°F) 

Define 

C SL - C C + XC AL (A ‘ 33) 

= the humid heat of the TEG-water solution, 

BTU/ (lb TEG) (°F) 

Differentiating equation (A. 32) 

dH L ‘ < C C +XC AL )dt L + C AL (t L ' V" 



(A. 34) 
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° r ’ dH L = C SL dt L = C AL <V t o )dX (A ' 35) 

Using equation (A. 12) 

(A. 36) 

L s 

Combining with (A. 31) gives the envelope II energy balance 

L S C SL dt I. * <h l a H <1Z - G S C AL dY)<t i- t L ) - G S H D' iY - Ua C (t L- t f >dZ < A ' 37 > 

Envelope III Energy Balance 

The overall enthalpy balance for control volume III is 
In Out 

L s\ + G S H G + Q D ■ L S (H L + dH L> + G S (H G + dH G> + G C (A. 38) 

Rearranging 

% m % - L S dH L ' °S dH G (A ’ 39) 

Substituting into equation (A. 37) 

L S C SL dt L = - G S IC s dt G + (C A < V t o ) - C AL (t L- t o )+X o +H D )dYl ' 

Ua c (t L -t f )dZ (A. 40) 

Interface Control Volume Energy Balance 
In Out 

W^WV'c’) + q SG VZ - * A V“( c *L (t r t o )> + q SL a H dZ (A - 41) 
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Rearranging, 

N A a M (C A (t G- t o > - C Al (t i- t o ) + V ' a H (h L <t r t L ) - h G (t G- t i» < A ' W > 

Solving for t^ 

t ± = WVVW.^VVlYc 1 (A. 43) 



Parameter Classification 

1. Known quantities for which inlet values are given 

v Y > v g > y* 

V x, H l , L, x, 

c f 



U, t 

o 

3. Variables with known functions 



V V V V V V y *> y m 

By starting at the top of the dehumidifier with the known quantities 
and solving the equations at each incremental step, dZ, the following 
quantities will vary; 

V Y - V G - y - c s’ 



"L ’ X ' V L ’ X - <W 
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V V V y *. V 



V V 

and will be obtained after each step and at the exit of the dehumidifier. 
Solution Sequence 

A solution sequence based on the differential section of Figure 
17 is: 



Y = 7^— 
1-y 

water mass ratio, 

X = - " 



1-x 



humid heat. 



liquid humid heat, 

C SL 

gas enthalpy, 

h g 

liquid enthalpy, 

=L + C SL ( V t o ) 



(A. 44) 



(A. 45) 



(A. 46) 



(A. 47) 



(A. 48) 



(A. 49) 



interface temperature, t^, from Eq. (A. 43); 

equilibrium vapor pressure above liquid, p*, from Eq. (2.2); 
mole fraction concentration in the gas in equilibrium with the bulk 
liquid from Henry’s Law, 




(A. 50) 
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mole fraction concentration in gas. 



'm 0.622 + 0.378 y 
mass transfer flux, 

\ ■ f g 



(A, 

(A, 



convection gas heat transfer coefficient corrected for simultaneous 
mass transfer 

h' = N A C A (A 

- 1-exp (-N A C A /h G ) 

This completes the supporting equations for the solution sequence, 
derivatives for the solution follow: 



dt G ‘ - h G a H ( V t l ) 



£ ■ d-«VH + c al ( - g s §»fvv + v- G s §>- 

dZ o oL 



(A 



(A 



dY N a 
dZ ^ ®M 
G s 

dX _ 1_. dY. 
dZ " L ^ S dz ; 



dy_ = 1 dY 

dZ= (1+y) 2 dZ 

1 



dx 
dZ : 



(A 

(A, 



(1 + X)" 



51) 

52) 

.53) 

The 

54) 

.55) 

.56) 

.57) 

.58) 

.59) 



dX 

dZ 



(. 60 ) 
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(A. 61) 




(A. 62) 




(A. 63) 




- 1 f > ■ "aN 



(A. 64) 



Before simulation of the contactor with the solution sequence can 
proceed, (1) input values must be initialized to correspond to the phy- 
sical system; (2) equations for fluid properties as functions of temper- 
ature and concentration are needed; (3) heat and mass transfer coeffici- 
ents for the physical system must be predicted or measured; (4) a ra- 
tional selection of the minimum number of variables necessary to specify 
the state of the model must be made. 

Equations (A. 54) through (A. 64) above are not all independent. 
Appendix B describes through bond graph techniques the selection of the 
minimum number of state variables to represent the LDD contactor. The 
results of Appendix B indicate that differential equationa (A. 54), (A. 55), 
(A. 56), (A. 63), (A. 64) above, plus necessary auxiliary equations, will 



model the contactor. 
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APPENDIX B 

BOND GRAPH CONTACTOR MODEL 

Introduction 

The bond graph method developed by Professor H.M. Paynter of MIT 
for systematic derivation of differential equations for the transport of 
power is described briefly: Bond graphs are a powerful analytical tool 

as they give a physical presentation of power flows by maintaining the 
topological structure of the system while allowing separate determin- 
ation of the sign of positive power flow and distinction between depen- 
dent and independent variables. While each engineering discipline has 
its own calculation methods, bond graphing is compatible with all dis- 
ciplines and is invaluable for model development at interfaces, i.e. 
electrical and fluid. Finally, the transformation from a bond graph to 
equivalent differential equations follows an exact method. 

Bond graph theory has been treated rather thoroughly elsewhere [13], 
[14], [25], [34], [35]. Thus, only those aspects necessary for under- 
standing the present model development will be discussed and only a 
meager sampling of the analytical potential of bond graph theory is pre- 
sented. 

General Bond Graph Model 

A bond graph representation of a control volume around the gas- 
liquid interface is shown in Figure 18. This bond graph is a distributed 
resistance-capacitance model of the interface and can be replicated n 
times to provide a model of any desired accuracy. Dynamic simulation in 
time and space is possible with this general model. 
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GAS, THERMAL 



LIQUID, THERMAL 




GAS, MASS 
TRANSFER 



LIQUID, MASS 
TRANSFER 




Figure 18. Bond Graph Representation of Coupled Thermal and Mass Transfer 
Operations in Liquid Desiccant Dehumidification 
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The zeros represent parallel junctions at which the sum of the 
heat or mass flows equal zero. A distinct temperature or chemical po- 
tential is associated with each zero junction. Two types of capaci- 
tances are represented on the graph: thermal, C and mass, C . There 

t m 

are two types of thermal resistance: convection resistance in the 

fluids, R t , and interfacial resistance, R^. Mass transfer resistance 

is also of two varieties: diffusion resistance, R,, and reaction re- 

d 

sistance, R . 

r 

The dotted line is a reaction coupling from the mass transfer por- 
tion of the bond graph to the heat transfer portion of the bond graph 
and represents the heat of condensation and dilution. 

Specific Bond Graph Model 

Looking at a time invarient solution of the general model allows 
the n replications to be done with a dynamic simulation routine if the 
differentials are taken with respect to dZ. Figure 19 shows this speci- 
fic simplified model with the addition of a cooling medium with heat 
capacity and temperature t^. 

Bond graphing allows a determination of the number of independent 
variables by counting the number of storage elements in the model. The 
temperatures t^, t^ and t^ are chosen to represent the mass transfer 
independent variables because of their physical significance and their 
association with the storage elements. Figure 20 is a word representation 
of Figure 19 and the contactor. (Note that the t's represent temperature 
differences. ) 
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At. 



'1 



At. 



'3 



At, 



cooling 

water 



air 




interface 




desiccant 



Q, 



■5 



*2 

heat source = fct(N A ) 



% 

heat source = fctCN^) 



air y_ 



interfacial x_ 

resistance N 
fct(t.) A 



liquid 



mass N a 



mass 



Figure 20. Word Bond Graph of Contactor 



Bond Graph Symbols 

On Figure 19 the bonds with arrows represent signal flows with no 
feedback along the same bond. For example, the interfacial resistance 
to mass transfer, R^, is dependent on the interface temperature. Solid 
bonds represent the power flow associated with heat or mass transfer. 

Ones represent series junctions at which the sum of the efforts, either 
temperature or chemical potential, is zero. 

In a true bond graph the product of the effort variable (above the 
bond) and the flow variable (below the bond) is power. Figure 19 is a 
pseudo bond graph because the concentrations are used as the mass trans- 
fer efforts rather than chemical potential and because heat flux is used 
as the thermal flow rather than entropy. With a pseudo bond graph the 
reaction couplings must be scrutinized as standard bond graph rules do 
not hold, but with consistent units, pseudo bond graphs follow the re- 
mainder of the rules for regular bond graphs. The pseudo bond graph was 
selected for the physical significance of the effort and flow variables. [14] 
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A "T" in a signal bond indicates the signal, t^ or N^, is trans- 
formed before use. The dotted line indicates the transformer modulus 
is controlled by t^. Linearity is not implied by the bond graph. Two 
mass flux controlled heat sources, and represent the latent heat 
of condensation and the heat of dilution. 

Heat flux in BTU/ (hr) (sq ft) and mass flux in lb water/ (hr) (sq ft) 
are the flow variables used in Figure 19. Table 1 lists the units of 
the thermal bond graph. 

Variable Units 

Effort °F 

Flow BTU 

(hr) (sq ft) 

Resistance (hr) (sq ft)(°F) 

BTU 

Capacitance BTU 

(hr) (ft) (°F) 

Table 1. Thermal Bond Graph Units 

Dimensionally, the product of resistance and capacitance is (ft) or 
the same as the solution increment. 

Differential Equation Model 

The differential equations which result from a rigorous application 
of bond graph techniques have been presented as equations (4.1) through 
(4.5) and are equivalent to equations derived in Appendix A. The com- 
puter simulation uses: (1) differential equations (4.1) through (4.5); 

(2) the auxiliary equation for t^, Eq. (4.6); (3) the equilibrium relation- 
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ship for vapor pressure over TEG, Eq. (2.2); and (4) the following se- 
quence of auxiliary equations: 
mass fraction in liquid, 




L 

percent weight fraction concentration of TEG, 



Henry's Law for the mole fraction concentration in the gas in equilibrium 
with the bulk liquid, 



P T = 100 (1 - x) 



(B.2) 



interface temperature from Eq. (4.6) 




equilibrium vapor pressure over TEG from Eq. (2.2) 



P* = fct(t i ,P T ) 




(B. 3) 



mass fraction in gas. 



G ~ G 



’S 



(B. 4) 



y = 



G 



mole fraction concentration in gas. 



y _ = z 

0.622 + 0. 378y 



(B.5) 



mass transfer flux (Fg is derived in Appendix C) 



1 - y* 




(B.6) 
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heat of dilution taken as 10% of the heat of condensation and sensible 
heat associated with the transfered mass, 

Q d = N a (100 + (t. - t 1 )) (B.7) 

mole fraction concentration in liquid. 



x = 
m 



x 

8.334 - 7.334 x 



(B.8) 



and humidity ratio. 




The computer model is now complete. 



(B .9) 
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APPENDIX C 

FLUID PROPERTIES AND TRANSFER PARAMETERS 



This appendix presents the essential but less than obvious relation- 
ships. The development of an equation for the vapor pressure over TEG 
from graphical data is discussed in Chapter II and is representative of 
the method of the numerous curve fits presented here. The remaining de- 
tails may be found in the computer program listing. Appendix E. 

Essential to the model is a numerical value for the mass transfer co- 
efficient, Fg. This coefficient is predicted from the more easily ob- 
tained gas convection heat transfer coefficient by means of the well 
accepted Chilton-Colbum analogy. Agreement of the analogy with newer 
data on mass transfer is very good [32], The mass and gas heat transfer 



coefficients are related by 



K p r 0.667 

F = — r - ) 
r n r v Co ' 



(C.l) 



Fg = mass transfer coefficient, (lb water) / (hr) (sq ft) 
where Cg = heat capacity of air at constant pressure, 

BTU/ (lb mass) 

Pr = Prandtl number 
Sc = Schmidt number 

h G = heat transfer coefficient between the gas and 
laminar TEG film, BTU/(hr) (sq ft) (°F) 

A numerical value for h^ is based on an empirical correlation for a 
representative finned tube extended surface heat exchanger reported by 



Rohsenow [30]. 
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0.134 k D G 

B . r max . 

n ' n ' 



0.681 0.333 

(Pr) (■ 



0.2 0.113 

) (?) 



where 

D = root diameter of tube, ft 
r 

k^ = thermal conductivity of air, BTU/ (hr) (f t) (°F) 



(C.2) 



G = mass flow rate at minimum cross section, 
max 

lb/ (hr) (sq ft) 

s = distance between adjacent fins, in 
1 = fin height or half the mean fin distance 
between tubes, in 
t = fin thickness, in 

y B = viscosity of air at bulk temperature, lb/ (hr) (ft) 



Eq. (C.2) shows excellent agreement with Kays [15] plate finned circu- 
lar tube data but which geometry differs from the modeled contactor. 

A curve fit developed for the Prandtl number for air is 



Pr = 0.721 exp(-196E-6 t g ) 

To compute the Schmidt number 

SC = ^~D~ 
v 

values for the viscosity, density and diffusivity are needed. 
Curve fits led to 



= 0. 0223(t ) 



492 [18.02 y + 28.96 (1-y)] 



(C. 3) 



(C.4) 



(C. 5) 
(C.6) 
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while the diffusivity is [2] 

0.000146 (t + 460) 2 ' 5 

n S- .... _ 

v t + 901 

g 



(C. 7) 



The thickness of a falling film of the water TEG solution on 
the vertical aluminum fin is [26] 




where 6 = film thickness, ft 

1^ = film liquid viscosity, lb/ (hr) (ft) 

T = liquid loading per unit width of fin, lb/(hr)(ft) 

Pp = density of film liquid, lb/(cu ft) 

2 

g = acceleration due to gravity, ft/hr 

The mass flow rate per unit width T can be calculated from the liquid 
recirculation rate but values for Up and pp for TEG-water solutions 
are obtained from curve fits 



1 4.46 E-7 (P_) 



[5.85 t 



(1.487 - 0.727- In P ) 



(C. 9) 



P D = 54.88 + 3.73 In P T - 0.0262 ^ (G.10) 

The liquid phase heat transfer coefficient was approximated as 
conduction dominated between the air-liquid interface and the liquid 



mid-point by 



h ^ 

\ = S 



(c.ii) 
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from the curve fit for the TEG-water solution conductivity 



k D = 0 



.789 - 0.14 In - 131E-6 t. 



(C.12) 



The heat transfer coefficient from the liquid mid-point to the 
aluminum fin was approximated by the empirical correlation for a 
falling film [17], [26] 

0.333 
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1.33 

\ P D C D 



L y, 



0.333 



M 



0.111 



(C. 13) 



where L = length of heat transfer surface, ft 
C D = C SL defined ln E< 1* ( A - 33 ) 

A curve fit for the specific heat of the TEG-water solution is 

= 0.4977 + 0.000444 t. + XC AT (C.14) 

D 1 AL 

A cautionary note with Eq. (C.13), [26] suggested a range of validity 
on L as 0.4 to 6 ft. Accordingly, the contactor height 1.12 to 2.24 
ft was used for L. This gave values for h^ from 40 to 60 while values 
for h^ ranged from 150 to 200. These values indicate that Eq. (C.13) 
may not be used since the heat transfer coefficient can not be less 
than that of pure conduction. Exclusive use of the conduction coefficient 
h^ is correct. The best approach is to calculate both h^ and h^ and use 
the larger for the entire liquid film. Sole use of h^ produces slightly 
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cooler operating lines as shown by Figure 21. In addition the typical 
overall heat transfer coefficient from the liquid to the cooling water 
increases from 14 to 20 while the typical overall heat transfer 
coefficient from the air to the cooling water increases from 6 to 7. 

The comparison in Figure 21 shows that the simulation results are not 
materially affected. 

A curve fit was developed from tube side water heat transfer co- 
efficient data [17] 

D - 2 6 

h t = (0.905 - 0.1863 In ( r - Q - - q - q^ ) ) (163 +1.731 t f ) V f * (C.15) 



where 

6 = tube wall thickness, ft 

w 

Vf = water velocity through tube, ft/sec 



The free flow area of Figure 8 is slightly reduced by the TEG-water 
solution film thickness so that the ratio of free flow area to face 
area is 



(C. 16) 

T u 

sh 

where P = number of fins per inch, 1/in 

T , = horizontal tube spacing, in 
sh 

The overall heat transfer coefficient from the gas to the cooling 
water is calculated by 
1 



U = J_ 
nh_ 



l 

nh T 



6 W , a H 
a k a^ h 



(C. 17) 



9100 
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Figure 21. Comparison of Conduction and Combination Liquid Phase 
Heat Transfer Coefficient 
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where 

2 3 

a^ = average tube surface area per unit volume, ft /ft 

= conductivity of tube wall, BTU/ (hr) (ft) (°F) 

5 = thickness of tube wall, ft 

w 

The calculations were done with the fin efficiency applied to both 
the fin area and the tube area. 

Applying the fin efficiency to the total area instead of just the 
fin area accounts for 97% of the heat transfer: 



r) _ . 66 

n a f a t ,66(.93)+.07 



(C.18) 



The remainder of the contactor geometry and transfer parameter 
calculations is straight forward. A complete listing may be found in 
Appendix E. 
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APPENDIX D 

TRI ETHYLENE GLYCOL PROPERTIES 

Triethylene glycol is a colorless, combustible, hygroscopic water- 
soluble liquid used as a chemical intermediate, solvent, bactericide, 
humectant, and fungicide. 

Physical Properties [12], [28] t 
Formula CH 2 - 0 - CH^ - CH^ - OH 
CH 2 - 0 - CH 2 - CH 2 - OH 



Molecular Weight 


150.2 


Boiling Point at 1 atm 


545.9 °F 


Vapor Pressure at 77°F 


<0.01 mm Hg 


Specific Gravity at 77°F 


1.119 g/cc 


Freezing point 


19 °F 


Pour Point 


-73 °F 


Viscosity at 77°F 


37.3 centipoises 


Surface Tension at 77°F 


45 dynes/ cm 


Specific Heat at 77°F 


0.53 BTU/(lb)(°F) 


Flash Point 


320 °F 


Fire Point 


330 °F 


Conductivity at 77°F 


0.13 BTU/ (hr) (ft) (°F) 
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Desired Desiccant Properties [33] 

A liquid desiccant selected for a dehumidif ication system should 

have the following characteristics: 

1. Suitable vapor pressure characteristics; should be capable of con- 
centration at low temperatures 

2. Noncorrosive 

3. Odorless 

4. Chemical stability; should not break down or vaporize in the range 
of interest 

5. Low viscosity and good heat transfer characteristics; must not 
crystallize or solidify near the operating range 

6. Nontoxic and nonflammable 

7. Widely available and inexpensive 

Triethylene Glycol Desiccant Properties [12], [28] 

1. TEG has a suitable vapor pressure curve and can be concentrated at 
temperatures as low as 120 °F 

2. TEG is compatible with steel if slight TEG discoloration due to iron 
contamination is acceptable. Otherwise stainless steel, aluminum, 
or resin coatings should be used. 

3. TEG is practically odorless. 

4. TEG is a superior choice considering both its hygroscopicity (item 1) 
and its high boiling point. 

5. TEG has a low viscosity and good heat transfer characteristics. 

6. TEG is low in chronic oral toxicity. It has been fed in the drinking 
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water of rats over a two year period at a 4% concentration with no 
apparent ill effect. Prolonged inhalation of saturated vapors of 
TEG has produced no ill effects on humans. Flash and fire points 
determined by the ASTM Cleveland Open Cup method for pure TEG are 
320 °F and 330 °F respectively. TEG presents a minimal fire hazard 
in storage or handling. 

7. TEG is a byproduct of automotive antifreeze. The cost of TEG com- 
pares to the cost of antifreeze. 
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APPENDIX E 



COMPUTER PROGRAM LISTING 
AND 



SAMPLE OUTPUT 



SEH = VORTHV R 511 21451 JOINT COMPUTER FACILITY 
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, X pL, * 2 



iH CNJ fO ITHD 1^- 00 O' 1—1 1- 



in Cl, H 
O U3 
W H C 



>-• PCI *C CM r- 



H [ i ?- H > 



►J *-} ,4 M M fvl ^ • 
O. «c «; >:: *. a: b-< t- 
SE.. (O U1 H H H «C h II 

Haiuioni^PHiS 



m m 



m 



n r» 

... • r- O CM 

(T • • (\ O H r 

II ii in ii • t -4 n 
a.: «« II tn li * sc 

X H a Cl< H & 111 

u O In Mm U 4 



ii ii m o cm 
U yi c, in m 
u« ^ r«j ii n 
EhOFmImO 



CM «- 

(N S. > 
00 iX (X, 

II CO CO 
1-3 D D 
E— ♦ E— • E— < 



WAPR0X = 4“. Y (18) = Y 

WK=8 . 1 
SH=. 01555 

GET CFM,DTA,FP,FPS,FT,GPMLIQ / LFN / PTEG / ROtfS,TF / TG,TL,TUSPH,TUoPV 



USER=VORTH*A 511 21451 JOINT COMPUTER FACILITY, MIT 
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t- E-i 
X 00 



X U .HI Cl. O GO i 
h fn h I/) II O i 



x cc r 

w c E 
H a: 
in * ; 
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O O 
3 X 03 
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: x to 



I d co H 
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O MI 
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3 D Q h M 3 



i- M m co i 
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: E-» M 
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Naa 
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* H E 
to X • 
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X X * 
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a a f 

E Ui c 
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II (0 II X II 



X • * «C ' 

x « x in ' 

O ?. I‘i »t H » 

E - 1 M « WOt 
* H O' H H E 

CO * oo + \ C 

3 h (O <«• E 

O O <»> X X * 

X El r- X X 



e E< Eh 



j II 



' ii : 



-> O II II to H < 
I X rf -< Eh X 
b co to x o *=c : 
-> Eh Eh X E-t X « 



IR0W=IFIX CROWS) 

FTU=LEN*TUT0T/12. 

FA = LEN * WI DTH /I 4 4 . 

END EXCHANGER DIMENSION CALCULATIONS 

BEGIN FLUID PROPERTIES CALCULATIONS 



USER=VOPTH«ft 511 21451 JOINT COYPUTFR FACILITY, ^IT 
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• to ii \ o 
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1 (S/FT) * * . 1 1 3/DI A 
E3=1 ./riCG 

DV=1.46F-4*(Y(2)+460.)**2.5/( Y(2) + 901 
SC=VGPS/( DGAS*DV ) 

FC = HCS * (?RG/5C)**.667/.24 
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